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Babu studied the neutrino spectrum obtained when one adds a charged singlet and a doubly charged singlet 
to the standard model particle spectrum. It was found that the neutrinos acquire a mass matrix at the two-loop 
level which contains one massless eigenstate. The mass matrix of Babu's model depends on an integral over the 
undetermined loop momenta. We present the exact calculation of this integral. 
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I. INTRODUCTION 

Zee observed that the addition of a charged singlet and addi- 
tional Higgs doublets to the Standard Model particle spectrum 
resulted in radiatively generated neutrino mass at the one-loop 
level Ql . It was consequently noted that if only one doublet 
couples to the leptons, the mass matrix takes a simple form 
and produces one mass eigenstate much lighter than the other 
two Ql- Babu studied the neutrino spectrum obtained when 
one retains Zee's charged singlet and adds a doubly charged 
singlet 1 3]. It was found that the neutrinos develop a mass ma- 
trix at the two-loop level, which contains one massless eigen- 
state, to lowest order. The mass matrix depends on an integral 
over the undetermined loop momenta of the two-loop diagram 
responsible for the neutrino mass, for which Babu gave an ap- 
proximate form. We present the exact calculation of this in- 
tegral. The analytic evaluation of related integrals, occurring 
when doublet neutrinos acquire mass at the two-loop level via 
the exchange of W bosons in the presence of singlet neutri- 
nos can be found in @]. 

In Sec. HTIwe briefly review Babu's model for the generation 
of neutrino mass. Sec. Mil contains the analytic evaluation of 
the integral which sets the scale for the massive neutrinos in 
Babu's model. The leading terms of this integral are presented 
in Sec. Hyjfor the relevant hierarchies of mass parameters 
involved. 



where f ab = -f ba and h ab = h ba . Gauge invariance pre- 
cludes the singlet Higgs fields from coupling to the quarks. 
The Higgs potential contains the terms: 

V(4>, h + ,k ++ ) = h i{h-h-k ++ + h+h + k—) + (2) 

which violate lepton number by two units and give rise to neu- 
trino Majorana mass contributions at the two-loop level (see 
Figure The neutrino masses are calculable and to lowest 
order the mass matrix takes the form: 

M ab = 8^/ ac ft, C(i m c m (i / C(i (/ t ) (ib , (3) 

where h ab = rjh ab with r\ = 1 for a = b and rj = 2 for a ^ b. 
m c .d are the charged lepton masses and: 
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cd ~ J Wfl (2^) 4 (p 2 - to 2 ) (p 2 - ml) (q 2 - to 2 ) 



(q 2 - to 2 ,) (p - q) 2 - ml 
If one defines: 

K cd = d, h Lh cd m c m d I cd , 

where no summation is implied by the repeated indices, the 
mass matrix may be written as M ab — {fKf') ab . Thus 
DetM = |Det /| 2 Detif = for an odd number of gener- 
ations (due to the anti-symmetry of /) and to lowest order the 



II. BABU'S MODEL 



Babu's model |3] includes two SU(2)l singlet Higgs 
fields; a singly charged field h + and a doubly charged field 
k ++ . The addition of these singlets gives rise to the Yukawa 
couplings: 
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FIG. 1: Two-loop diagram responsible for neutrino mass in the Babu 
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neutrino spectrum contains one massless state. In what fol- 
lows, the exact calculation of the integral I c d is performed and 
the asymptotic behaviour for the cases m k ^> m/j 3> m c , rrid 
and mh ^S> m k 3> m c , m d is presented. Note that the hier- 
archy m/j,TOfc ^> m c ,m<j is a phenomenological constraint, 
whilst the relative size of m k and mt is not predetermined. 



III. EVALUATION OF I cd 

After performing a Wick rotation and letting q — > — g, the 
integral / c( j may be written as: 
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(q 2 + m 2 d ) (p + q) 2 + ml 

where the momenta p and q are now Euclidean four-vectors. 
We adopt the notation of |7] and define: 

(M 11 ,M 12 ...M lni \M 21 ,M 22 ...m 2n2 \M 31 ,M 32 ...M 3na ) 
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so that 7 cd = (5^p-("i/i, m c \m h , m d |m fc ) when n = 4, where 
7i is the space-time dimensionality. Note that all momenta 
in are Euclidean whilst the definition used in |7] contains 
Minkowski vectors. It is possible to express I c d as a linear 
combination of integrals with less than five propagators. One 
may use partial fractions to obtain relations like: 
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{(to |toi|to 2 ) - {m\mi\m 2 )} . 



At this point one may proceed to evaluate the integral 
(mo | To! | rri2 ) to determine I c d but it is preferable to express 
I c d as a combination of integrals of the form (2too|tt7i \m 2 ) 
(where (2to-o|toi \m 2 ) is shorthand for (too,TOo|toi|to 2 )). If 
one evaluates integrals of the form (too|t77i|to2) by intro- 
ducing Feynman parameters, some of the ultraviolet diver- 
gences are transferred from the radial integrals in momentum 
space to the Feynman parameter integrals. It is easier to han- 
dle the divergences when evaluating (2mo|mi|m 2 ), as they 
may be completely contained in the momentum space inte- 
grals. Integrals of the type (too|toi |to 2 ) have been consid- 
ered in M l a] an d single integral representations have been 
obtainedlUMlllil. 

To express (too|toi|to2) in terms of integrals of the 
form (2to(i Itoi \m 2 ) one may use the partial p operation of 
't Hooft 111 311 . Essentially one inserts the identity expression: 



1 " 

— y 

2n ^ 



i=l 



dpi 
dpi 



(8) 



into the integrand of (too|toi |to 2 ), performs integration by 
parts and rearranges the resulting expressions to obtain the 
relationship: 



(to |toi|777 2 ) = x {too(2to |toi|to 2 ) 

3 — n 

+TO 2 (2777 1 |too|?772) + 7772(27772177701777!)} . 



(9) 



which may be applied to I c d to give: 
1 1 1 



{{m c \m d \m k ) 



(2tt) 8 (to 2 - to 2 ) (to 2 - to 2 ) 
(m h \m d \m k ) ~ (m a \m h \m k ) + {m h \m h \m k )} . (7) 



We note that equation (|SJl is defined for integral 77, but the 
resulting relationship may be analytically continued to non- 
integral 77. Using equation (|9jl on I c d gives: 
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x { to 2 [(2m c \m h \m k ) - (2m c \m d \m k )} + to 2 , [(2m d \m h \m k ) - (2TO d |m c |m fc )] 



to 2 [(2TO fc |TO/ l |?77 £ /) + (2TO fe | to c \m h ) - (2m k \m c \m d ) - (2m k \m h \m h )\ 



+ ml [(2m h \m d \m k ) + (2m h \m c \m k ) - 2(2m h \m h \m k )] } 
Thus evaluation of the generic integral (2to|toi|to2) allows one to determine I c d- We obtain: 

7T 4 (7TTO 2 )"- 4 r(2- in) f 1 
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where 



where the function f(a, b) is given by: 



ax + 6(1 — x) 
x(l — x) 



rn 
m 



This result is in agreement with the result for Q(m,m±, m 2 ; 0) 
in pill and differs by an overall minus sign to that obtained 
in |3I due to our different definition of (2to|toi |ma) in terms 
of Euclidean momenta. Letting n = 4 + e and expanding for 
small e gives: 

(2m|mi|m 2 ) 

^ + -(l-2 7B -21og(^m 2 )) 

1 1 

~2 ~ 12 

.2/ 2 



2 ~ ~ 7 s + (1 - 27e) log(^m 2 ) 

log 2 (7Tm 2 )-.f(a,b)]+0(e) (11) 
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and the dilogarithm function is defined as: 

log(l - y) 



Li 2 (x) 



Evaluation of /(a, 6) gives: 
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We have introduced: 



x\ = -(1 + t- aH 
yi = -(1 + a- 6- 



^2 = ^(1 
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and: 



(l-2( a + 6) + (a-6) 2 ) 1 / 2 , 



where under a ^> 6 we have ^> 7/j for z = 1, 2. The expres- 
sion dl2l can be shown to be equivalent to: 



f(a,b) = -ilogalog& 
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(13) 



which is the symmetrized version obtained in 0]. For nu- 
merical evaluation via mathematica etc, it may be more con- 
venient to use the form dl2l . If one uses Jl 31 for arbitrary 
masses it is possible for the logarithms to go negative, thus 
producing a non-zero imaginary component for f(a,b). By 
judicious logarithmic branch choice one can ensure that the 
imaginary components of f(a, b) cancel, whilst the real part 
of /(a, b) is independent of the branch choices. This real part 



is always in agreement with the form dl 2i of f(a, 6), but in 
the latter the imaginary parts cancel when using the principal 
branch for all dilogarithms. The result for I cc [ will ultimately 
be a combination of functions of the form f(a,b). As I c d 
is evaluated in the neutrino rest frame, it is not possible to 
cut the two-loop diagram it represents and obtain an energy- 
and momentum-conserving sub-graph (ie the neutrino at rest 
doesn't have enough energy for the internal particles to be real 
and on the mass shell). Thus the amplitude for the graph must 
be real 1 15], forcing the function f(a, b) to be real. Equation 
(II It together with either dl 31 or fl!2i gives the final result for 
(2m|mi|m 2 ). 

I c d is now obtained by substituting the result for 
{2m\mx\m 2 ) into the expression dlO> . Inspection of dl 01 
shows that we may use (2m|mi|m 2 ) — > 7r 4 /(a, b) when eval- 
uating I c d, as the constants and logarithms with massive ar- 
guments occurring in (2m|mi|m 2 ) cancel amongst the terms 
with a given mass coefficient. 



IV. DOMINANT BEHAVIOUR OF I cd 

We now find the asymptotic behaviour for the cases 
nik ^> nih and m/i ^> m^. The dominant terms are found by 
using the expansions for (2m|mi|m 2 ) presented in |7]. Utilis- 
ing the symmetry properties of (2m|mi |m 2 ) and the structure 
of the expression (II 01 for l c d allows one to obtain the leading 
terms relatively easily. 
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1. m k > m h . 



m h » m k . 



In the expression for I c d, equation (I10> . we may safely 
neglect the terms with lepton mass coefficients. The terms 
with coefficient m 2 have the form (2m k \m x \m y ) and the 
expansion for (2m\M\\M2), where m ^> Mi,M 2 , given 
in Q, may be used. All constants and logarithms of m 2 
may be ignored as they cancel amongst the various terms. 
The expansion is given in terms of a = (A/i/m) 2 and 
b = (M 2 /m) 2 . By noting the form of (TTOl and that 
(2m | Mi | Mi) = (2m | M 2 1 Mi), it is seen that in the expan- 
sion of a term (2m k \m x \m y ), any term depending only on a 
or only on b will occur in the expansion of another term with 
a relative minus sign. So one need only retain terms contain- 
ing both a and b. Any terms which contain a lepton mass 
will be suppressed, meaning the term (2m k \mi l \irih) dom- 
inates. Defining h = (mh/m k ) 2 and using (2rn\Mi\M 2 ) 
in 1 7], keeping only the terms containing both a and b, gives: 

(2m k \m h \m h ) = -tt 4 {h 2 + 6h 3 ~ 2 \og(h) {h 2 + 5/i 3 } 



-[log 2 W + y]{/ i 2 + 4/ l 3 } + ....j 



(14) 



All the terms with a coefficient m 2 , in I c d have the form 
(2mh\m x \mk) and will be expanded in terms of h and 
(m x I m k ) 2 ■ The form of dlOl shows that any term in an ex- 
pansion depending only on h (not including the cross terms 
when m x — mh) will occur in another expansion with a rela- 
tive minus sign. If m x is a lepton mass, the leading contribu- 
tions from (2m,h\m x \mk) will contain factors (m x /m k ) 2 and 
consequently be suppressed. Thus to leading order: 

m 2 {(2m h \m d \m k ) + (2m h \m c \m k ) - 2(2m h \m h \m k )} 



-2m 2 h (2m h \m h \m k ) 



(15) 



Using the expansion of (2Mi\M2\m) for m ^> M\,M2 
from [7] and noting the above gives: 

(2m h \m h \m k ) = -vr 4 - ^-h 2 + y (h + 2h 2 ) 

+ log 2 (/i) {h + 3h 2 }+ log(h) j/i + y/i 2 j + ..1(16) 

as the leading terms which contribute to I c d- Combining the 
expressions for (2m/ l |m; l |mfe) and (2m k \irih\mh) and noting 
the appropriate coefficients from equation i fTOl gives: 
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where the lepton masses have been neglected relative to the 
scalar masses in the factors (m, 2 — m 2 ) _1 (m 2 — m 2 ,) -1 . 
Babu's approximate form reproduces the log 2 term of d!7l > in 
the h — > limit. It is necessary to retain both (2m/ l |m? ! ,|mfc) 
and ^niklmhlmh) to obtain the leading terms of |17t . The 
leading terms of I c d may be obtained by setting m c = = 
but one may not take m,h = when m k ^> rrih as I c d devel- 
ops an infra-red singularity. This singularity is seen to mani- 
fest itself in the leading terms of equation < l 1 71 as a logarithmic 
singularity in the limit m/, — > 0. 



The terms with lepton mass coefficients may again be ne- 
glected to obtain an expansion in terms of k = {rn 2 ,/rnf l ). The 
presence of the factors (m 2 — m 2 ) _1 (m 2 1 — m 2 ,) -1 w m^ 4 
in fl!0l > means that the terms with coefficient m 2 will be sup- 
pressed relative to those with coefficient m 2 . The expan- 
sion of (2m|m|M) for m > M in 0] may be used for 
(2m? l |m/i|m/ c ) where only terms originating in the f(a,b) 
portion, as defined in dl 11 1 need be retained. The result is: 



(2m h \m h \m k ) 



log(fc) 



36 
1 , 



12 



-k 1 



■(18) 



For mt ^> m k , contributions from the terms (2rrih\m x \m k ), 
where m x is a lepton mass, must be included. The lead- 
ing terms are those which do not contain factors (m 2 /m?_). 
The expansion for (2m\M 1 \M 2 ) with m > M X ,M 2 in jj 
is again employed. Retaining only those contributions which 
stem from f(a,b) gives: 



(2m h \m x \m k ) 
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with this expression holding for m x = m c and m x = rrid. 
Using dl 81 and fl!9l > and neglecting terms with coefficients 



m" c , m 2 and m 2 in i ll Oi l, gives the result 
11 
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Thus for large mi the leading term is: 



(4tt) 4 - 
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O -r . (20) 



(21) 



The m h factor is expected on dimensional grounds and in 
this limit the coefficient of m^ 2 contains no logarithmic sin- 
gularities. We see that the dominant terms for I c d, when 
rrih ^S> m k , are found by considering the terms with coef- 
ficient mi in JlOb . At higher order in (m k /mh) 2 contribu- 
tions from the terms with coefficient m\ will also become im- 
portant. Note that \2\\ is independent of m k and the lepton 
masses and that its form may be understood as follows. De- 
fine: 



h{k 2 ) 
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(p 2 — m?) (p — fc) 2 



(<z 2 



2) {k - q) 2 {p _ q) 2 



(22) 



This corresponds, up to multiplicative factors, to I c d with 
rn c = rrid = m k = and m/j — > m when the external mo- 
menta k is taken to be zero: 

k - {2 lr /2(fc2) 



I cd\ mc=md=mk=0 - lim 
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The integral hik 2 ) has been evaluated 
suit: 



1 



], leading to the re- 



■M me= „i d=nH= o - ( 47r )4 m 2 C( 2 )- ( 23 ) 

Noting that £(2) = ^ and comparing A23I with fl21i it is 
seen that the leading term of I c d, for rrih ^> m^, is simply the 
result one obtains for I c d when the lepton masses and rrik are 
taken to be zero (observe that for rrih ^> we may safely 
set m c , mrf and to zero to obtain the leading term as I c d 
remains infra-red finite in this limit). 

After completing this work we became aware that the ex- 
plicit asymptotic terms in equations dl7l and M\\ were given 
in 07 



V. CONCLUSION 



The analytic result for the integral: 

d 4 p f d*q 1 
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(2^) 4 J (2^) 4 (p 2 - m 2 h ) (p 2 - m 2 ) (q 2 - m\) 
1 1 

(24) 



(q 2 - m 2 d ) (p - q) 2 - m\ ' 



which appears in the lowest order neutrino mass matrix in 
Babu's model, has been obtained. The result comprises of 



equation dlOi together with ( II It and either Jl 31 or dl 2i for the 

terms (2rri\rri\\rri2) . We note that for rrik 3> rrih the leading 
term is: 



led 



_J L log 2 ( 

(^ m 2 k l0g \^ 2 



The term in the asymptotic expansion was extracted and 
is in agreement with the result of 1 16]. The leading term when 
rrik 3> rrih is independent of the lepton masses but requires 
a non-zero value of rrih to avoid an infra-red singularity. For 
the reversed hierarchy, rrih ^ rrik, the leading term is: 



1 



2 TT 2 



(4tt) 4 to 2 6 



which is the result one obtains for I c( i upon setting the lepton 
masses and to zero and is in agreement with that obtained 
in ill- 
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